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ABSTRACT 

We present ultradeep optical spectroscopy obtained with FORS2 on VLT of seven Lyman-break 
galaxy (LBG) candidates at z > 6.5 selected in the GOODS-S field from Hawk-I/VLT and 
WFC3/HST imaging. For one galaxy we detect a low significance emission line (S/N < 7), lo- 
cated at A = 9691.5 ± 0.5 A and with flux 3.4 x 10 _18 erg cm -2 s -1 . If identified as Lya, it places the 
LBG at redshift z = 6.972 ± 0.002, with a rest-frame equivalent width EW rf = 13k. Using Monte 
Carlo simulations and conservative EW distribution functions at 2 < z < 6, we estimate that the 
probability of observing no galaxies in our data with S/N > 10 is ~ 2%, and that of observing only 
one galaxy out of seven with S/N = 5 is ~ 4%, but these can be as small as ~ 10 -3 , depending on 
the details of the EW distribution. We conclude that either a significant fraction of the candidates 
is not at high redshift or that some physical mechanism quenches the Lya emission emerging from 
the galaxies at z > 6.5, abruptly reversing the trend of the increasing fraction of strong emitters with 
increasing redshift observed up to z ~ 6.5. We discuss the possibility that an increasingly neutral 
intergalactic medium is responsible for such quenching. 

Subject headings: galaxies: distances and redshifts - galaxies: high-redshift - galaxies: formation 



1. INTRODUCTION 

During the past year, a suite of new near-infrared 
(NIR) surveys has extended the search for star-forming 
galaxies to redshift 6.5 < z < 10 using t he well-proven 
Lyman-break technique ( Giavalisco|2002| and references 
therein). With respect to lower redshift, the number 
density of UV-selected galaxies decreases (e.g 
et al.||2009[ IMcLure et al.||20l0l [Castellano et al 



Ouchi 



Bouwens et al.||2010a| |WiIkins et al. | |2010bp , their U 



energy distribution (SED), surface brightness and mor- 
phology, with the result that the measure of the luminos- 
ity function remains subject to uncontrolled systematic 
errors. 

In practice, given the marked decrease in sensitivity of 
current spectroscopic observations at increasing redshift, 
the spectral confirmation of ga laxies at z > 5 relies heav- 
ily on their Lya emission line, ( [Stark et al.|2010||Vanzella| 
et al.||2009| S10 and V09 in the following). Indeed, red- 



continuum becomes bluer implying either reduced dust 
obscuration or poorer metal enrichment or both (e.g. 



shifts derived without Lya typically have lower confi 



Finke lstein et al.|2009| |Bouwens et al.|2010b| |Schaerer fc| 



de Barros|2010 1, and their stellar masses are, on average, 
small er than those of th eir lower redshift counterparts 
(e.g. " 



dence, although their number may be comparable ( Dou 
glas et al.|2010| D10 in the following). The line in itself is 



an important dia gnostic of physical processes at work in 
the galaxies (e.g., G iavalis co et al.|1996||Pentericci et al 



Labbe et al.||2010p. Unfortunately, these results |2010| |Shapley et al.| |2003p, since its strength and veloc 



are based on color-selected samples with no spectroscopic 
validation. At the time of writing, spectroscopic detec- 
tions of onl y a few individual objects haye been obtained 
at z > 6.6 dlye et al.|2006l|Greiner et al. | 2009 [ |Salvaterra 



et al.|2009|]Tanvir et al.||2009) 



'1'he lack of knowledge of the true redshifts of the cur- 
rent z ~ 7 candidates places significant limitations on our 
ability to robustly measure the properties of the galaxies 
at this critical cosmic epoch. For example, the fraction 
of interlopers and the redshift distribution of the sam- 
ple galaxies are necessar y to robustly measure t he UV 
luminosity function (e.g., Reddy & Steidel 2009). Cur- 
rently, the former remains unknown, and the latter is 
estimated with Monte Carlo simulations under various 
assumptions for the intrinsic distributions of UV spectral 
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ity profile depend on the instantaneous star-formation 
rate, dust content, metallicity, kinematics and geometry 
of the interstellar medium. Particularly relevant here is 
the evidence that the fraction of Lya emitters i n UV- 
selected sample s increases with redshift (V09,S10 |Reddy| 
fc Steidel|2009}|Stanway et al.|2007| S07 in the following) 



and that the traction of galaxies with a large Lya equiv- 
alent width (EW) is substantially larger at fainter UV 
luminosities. 

Finally, the very visibility of the Lya line during the 
ending phases of the cosmic re-ionization is subject to 
the damping effect of a n increasing neutral intergalac 



tic m edium (IGM) (e.g., Zheng et al. 2010 Dayal et al 



20101, expected to attenuate most ot its luminosity and 



make the earliest galaxies consequently more difficult to 
identify. Hence, the line profile and the evolution of its 
EW are sensitive diagnostics of the ionization state of 
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the high redshift IGM. 

To address these issues we have started a campaign 
of spectroscopic follow-up of z ~ 7 "Z-dropout" can- 
didates, selected from high-quality imaging surveys ob- 
tained with VLT/Hawk-I and HST/WFC3. In this paper 
we present the fir st results from a sampl e selected in the 
GOODS-S field ( |Castellano elTaT1|2l)iol CIO in the fol- 
lowing). All magnitudes are in the AB system, and we 
adopt H = 70 km/s/Mpc, ti M = 0.3 and fi A = 0.7. 

2. TARGETS AND SPECTROSCOPIC DATA 

This initial spectroscopic sample includes relatively 
bright Lyman break galaxy (LBG) candidates at z > 6.5 
(listed in Table ]!]), five from the H awk-I images (4 from 
CIO and 1 from|Hickey etaT1(|2010|)) and two from WFC3 



( Oesch et al.||20To[ |Wilkms et ai.||2010a[ ), spanning the 
magnitude range Y ~ 25.5— 27.5. We tilled empty slitlets 
in the multi-object slit masks with other candidates of 
lower quality and / or at lower redshi f t, incl uding a candi- 



date brown dwarf (Mannucci et al. (20071, C10) and i— 



dropouts from the GOODS survey not observed by V09 

2.1. Observations 

Observations were taken in service mode with the 
FORS2 spectrograph on the ESO Very Large Telescope, 
between 12 November 2009 and 14 January 2010. We 
used the 600Z holographic grating, that provides the 
highest sensitivity in the range 8000 — 10000A with a 
spectral resolution R ~ 1390 and a sampling of 1.6A per 
pixel for a 1" slit. The data presented here come from 
the coaddition of 75 spectra of 842 seconds of integration 
each, on a single mask, for a grand total of 63150 s (17.5 
hr), with median seeing around 0.8". Each slitlet was 
1" wide and 14" long, to maximize the number of slits 
available while allowing a careful sky subtraction. All 
our high priority targets were placed at the center of the 
slits, and spectra were taken in series of three different 
positions, offset by ±2" in the direction perpendicular to 
the dispersion. 

Since our objects are extremely faint, the slit centering 
was based on the astrometry solution obtained from the 
Hawk-I images, which is which is well aligned to the ACS 
one. We have directly verified this by placing a few bright 
objects from the ACS catalogs in small slits, and ensuring 
that they were correctly aligned during the observations. 
It is also reassuring to note that three faint i-dropouts se- 
lected from the ACS catalog which were placed in slitlets 
using the same astrometry, have a clear Lya detection at 
z ~ 5.94 (full details will be given in a future paper). 

Data reduction was performed using an optimized ver- 
sion of the recipes adopted in V09 and previous papers. 
After standard flat-fielding and wavelength calibration, 
we subtracted the sky emission lines with two different 
procedures. In the first case (Polyn in the following) we 
fit polynomials of order n (from 1 to 5) to the sky inten- 
sity at each pixel position. This procedure in principle 
ensures the highest S/N, but is sensitive to systematics 
induced by defects in the detector or in the slit manufac- 
turing. A safer but somewhat noisier approach (ABBA 
in the following) is to subtract the sky background be- 
tween two consecutive exposures, exploiting the fact that 
the target spectrum is offset due to dithering. We found 
the spectra obtained with the two techniques entirely 



TABLE 1 

-DROP TARGETS IN GOODS-S 



ID 


R.A. (deg) 


DEC. (deg) 


Y 


Z-Y 


M (a) 
M uv 


G2_1408 i > 


53.177382 


-27.782416 


26.37 


>2.1 


-20.49 


G2_2370 b 


53.094421 


-27.716847 


25.56 


1.68 


-21.27 


G2_4034 i > 


53.150019 


-27.744914 


26.35 


>2.1 


-20.50 


G2_6173 i > 


53.123074 


-27.701256 


26.53 


>1.9 


-20.33 


H_9136 c 


53.072574 


-27.728610 


25.90 


1.29 


-20.94 


W_6 d 


53.100392 


-27.703847 


26.93 


1.17 


-20.38 


0_5 e 


53.177376 


-27.7920551 


27.52 


1.61 


-19.67 



a - Computed at z = 6.8 

b -I Castellano et al. (2010), Yqpen Hawk-I 
c - [Mickey et al |(|2010[, Y qpf.n Hawk-I 
d - jWilkins et^t[l|2010a| , Y 98 WFC3 - ERS 
e - [Oesch et affpOlo} , Y w5 WFC3 - HUDF 
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Fig. 1. — Lower: Limiting flux density (at ler level) resulting 
from our observations. Upper: Corresponding 10cr limit on the 
rest-frame equivalent width of a Lya emission line as a function of 
redshift. Colors and line widths correspond to different observed 
magnitudes in the Y band, as shown in the legend. 

consistent. Finally, spectra were flux-calibrated using 
the observations of spectrophotometric standards. Slit 
losses are small, given the extremely compact size of the 
targets and we neglect them in the subsequent analysis. 

The r.m.s. of the resulting spectra, which will be used 
later to determine the probability of our results with a 
Monte Carlo simulation, has been estimated "by first 
principles", i.e., computing the absolute r.m.s. of each 
frame from its raw counts c as yjcfg) (where g is the 
e~-ADU conversion factor) and propagating it through 
all the reduction steps. It turned to be in excellent agree- 
ment with the observed r.m.s. in the region between 
8150 and 8250A that is devoid of sky lines and with 
the predicted efficiencies estimated by the ESO exposure 
time calculator. The resulting la limiting flux density is 
shown in the lower panel of Fig[T] 

To obtain the corresponding limit on the detectable 
EW for a Lya line, we have computed three differ- 
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Fig. 2. — Spectrum of the candidate G2_1408, showing a tenta- 
tive emission line at 9691. 5A. The two upper panels show the 2-D 
spectrum of the sky emission and of the sky-subtracted object, as 
indicated. The x-axis is in wavelength, in the same range of the 
three spectra below. The 2-D spectrum of the galaxy has been 
divided by the r.m.s. to remove obvious spikes due to bright sky 
lines, and slightly filtered with an adaptive mesh. The three 1-D 
spectra in the lower part show the extracted spectrum (over 4 pix- 
els) with the two different techniques for sky subtraction, and the 
sky emission at the same wavelengths (see legend) . In these panels 
the spectra have not been divided by the r.m.s., nor filtered. 

ent cases, assuming continuum magnitudes of m = 
25.5,26.5,27.5, to span the luminosity range of our tar- 
gets. For the computation we assume that the flux profile 
is a Gaussian with FWHM= 10A. The resulting limiting 
EW is shown in the upper panel of Figjl] computed at the 
lOcr level. We could detect weak [EW ~ 5A) Lya lines 
in our brightest galaxies, and even for the faintest ones 
we are able to reach EW ~ 50A over a significant frac- 
tion of the redshift interval. This range of sensitivity is 
similar to that of z ~ 5-6 surveys (S07, V09, D10, S10). 

2.2. Results 

We detect only one weak emission line, centered at 
9691.5 ± 0.5A in the spectrum of the object G2_1408. 
This galaxy is the brightest candidate identified in the 
Hubble Ultradeep Field (HUDF) area, and o ne of the 
bright e st in CIO. It was first detected by |Bouwensl 



et al. (2004) in the NICMOS HUDF data, and subse- 
quent lflclehtified also by C10 and in the HUDF WFC3 
data dBouwens et al. 2010a lOesch et al. 2010; McLure 



'5 



et al. 2010| |Bunker etTaTf 2009" U From the clear elon- 
gation observed in the WFC3 images, one can exclude 
the possibility that it is a brown dwarf. The 2-D and 
1-D spectra of G2_1408 are shown in Fig(2] The spec- 
tral feature is extended over 4 pixels in the spatial di- 
rection, consistent with the average seeing. The FWHM 
is ~ 10A, significantly larger than any feature due to 
noise. The weak emission line has a total observed flux 
of 3.4 x 10" 18 erg cm. _2 s _1 . The formal S/N is 7, but this 
estimate does not include systematic errors, and should 
be considered as an upper limit. We made extensive tests 
to verify the reliability of this detection. We verified that 



the feature is present both in the Polyn and in the ABBA 
reductions, as shown in Figj2j We then inspected all the 
75 individual spectra to ensure that the feature is not 
due to an artifact, and that it is still detected when we 
separately summed the data in two halves. Because of 
the large color break (z — Y > 2.1) measured in the 
HUDF data and the non-detection in the BVI bands, 
an identification of this line with a lower redshift [Oil] 
or Ha would imply a very peculiar SED, unlike that of 
currently known galaxies. This cannot be excluded a 
priori. 

We note that there is no evidence of the asymmetry 
that is expected (but not required, see discussion below) 
for a z ~ 7 galaxy, although the S/N is too poor to reach 
any firm conclusion about this. 

Based on these tests, we conclude that the feature is 
likely real and due to Lya emission from a z — 6.972 
galaxy (z — 6.970 if computed at the blue edge of the 
line), although this should be validated by independent 
and possibly deeper observations. No continuum is de- 
tected in the spectrum: if we estimate it from the Hawk-I 
Y-band magnitude (Table 1), the line flux translates into 
an observed EW of 103A, corresponding to 13A if placed 
at z = 6.972. 

We do not identify any other emission lines from ob- 
jects in our sample. We only detect a faint continuum 
from two objects, namely G2_2370 (the brigh test in our 
sampl e ) and the brown dwarf candidate of Mannucci] 
et al. ( |2007 ). In both cases, the continuum is consis- 
tent with the broad-band magnitudes but the low S/N 
prevents us from deriving any robust information about 
their spectral type or redshift. 

3. THE EXPECTED NUMBER OF Lya DETECTIONS 

The key result of our observations is the lack of promi- 
nent Lya emission lines in our sample, which may im- 
ply a rapid evolution in the physical properties of z > 6 
galaxies and/or in the surrounding IGM. To quantify this 
issue, we have carried out the following Monte Carlo sim- 
ulations under the assumptions that a) all our 7 candi- 
dates are indeed z ~ 7 galaxies; and b) the distribution 
of the Lya intensity in galaxies as a function of their 
rest-frame continuum magnitude Mjjy does not change 
significantly from z = A — 6 to z = 7. 

For the redshift distribution expected for our sample 
we use the result by C10 (see their Fig 7), which has 
a broad maximum from z = 6.4 to z — 7.1 and tails 
that extend to z = 6 and z = 7.5. The distribution of 
the Lya intensity in galaxies at z = 3 — 6 has been in- 
vestigated in a number of studies (S07, V09, S10, D10), 
showing that the intensity of Lya is anti-correlated with 
rest-frame UV luminosity. No measure of the depen- 
dence of the EW distribution as a function of Muv nas 
been obtained, however. We model the EW distribution 
assuming that at EW> it is represented by a Gaus- 
sian centered on EW=0 with an additional constant tail 
up to 150A, and at EW< by a constant level down to 
some EW m i„ value, and null below. We take the width 
of the Gaussian and the two tails to reproduce the re- 
sults of V09 and S10 at different rest-frame magnitudes. 
Specifically, we derive from the bright galaxies in V09 a 
standard deviation for the Gaussian of 10A, and assume 
that it is constant at all magnitudes. We then divide 
our sample in two luminosity bins (—20.5 < Mjjv and 
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—20.5 < Mjjv < —19.5) and adjust the two tails in or- 
der to reproduce the fraction of galaxies with EW > 50A 
given by S10 and the fraction of galaxies with EW > 5 
and EW> 20 A (for the two bins, respectively), as given 
by the V09 data. The resulting distributions are shown 
in Figure [3] for the two magnitude bins, and are reason- 
ably similar in shape to the EW distribution at z ~ 5 — 6 
(S07), and show a m oderate evolution from the z ~ 3 — 5 
( jShapley et alp003| D10) one. 

We then compute the probability of detecting N 
Lya lines at a given S/N in our sample of 7 objects. For 
each object we randomly extract a redshift from the CIO 
distribution, we compute the corresponding Mjjy from 
the observed Y band magnitude (taking into account the 
IGM absorption at that redshift), and we then randomly 
extract an EW from the corresponding distribution. If 
the EW is larger than the minimum detectable EW at 
the corresponding wavelength (Figfll) for a given S/N we 
conclude that the object would be aetected. W e assume 
FWHM=10A for the line, as found at z = 6.9 by |lye et al." 
(2006) (see also Fig [I]). Clearly, intrinsically broader lines 
would be harder to detect. We perform this exercise 10 5 
times over the whole sample, requiring S/N> 10 for the 
detection (larger than the S /N of the possible detection 
in G2_1408), and we finally obtain the probability distri- 
bution shown in the lower panel of Fig. [3] Under these 
assumptions, the probability of detecting no Lya line in 
our sample is very small, about 2%, while the typical 
number of Lya that we should have detected is between 
2 and 4. We also find a low probability (4%) of having 
1 detection at S/N> 5, as found in our sample. The 
same probability adopting the S07 distribution would be 
much smaller (~ 10~ 3 ), because of the s ubstantial tail 
of obj ects with large EW. Even using the Shaplcy ct al. 
(2003) distribution, which has a lower fraction of high 
EW objects, the probability is still rather low (9%). We 
conclude that, with all the obvious caveats due to the 
small size of our sample and to possible observational 
mishaps, the lack of prominent Lya lines in our sample 
is statistically significant. 

4. IMPLICATIONS AND DISCUSSION 

On a practical level, our results show how challeng- 
ing it is to obtain large samples of spectroscopically con- 
firmed galaxies at z > 6.5 with current instrumentation, 
especially if one aims at reaching the level of complete- 
ness 50%) needed to robustly measure the luminosity 
function. Our observations imply that this goal will have 
to wait until a future generation of instruments is avail- 
able, either 8m telescopes equipped with multi-object 
spectrographs more efficient in the z and Y bands or, 
more likely, the new generation of telescopes, such as the 
James Webb Space Telescope or 20-40m ground-based fa- 
cilities. Nonetheless, our analysis appears to show that 
the failure to detect prominent Lya in our sample is not 
only due to the insufficiency of current instrumentation. 

One possibility is that a significant fraction of the can- 
didates are lower redshift interlopers. We test this pos- 
sibility by extrapolating to z ~ 7 the observed contami- 
nation in spectroscopic samples at z ~ 4, 5 and 6 (V09, 
Table 4 of B, V and i dropouts), which increases with 
redshift. We assume that amongst our z-band dropout 
sample the fraction of contaminants could be ~25%, i.e., 
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Fig. 3. — Simulations on the expected number of Lya emitters in 
our sample. Upper: The adopted distribution of rest-frame EW for 
the two extreme values of rest-frame UV luminosities in our sam- 
ple (red a nd black continuous histogram) . The shaded histogram 
shows the |Shapley et al.| l |2003| distribution at z ^ 3 and the blue 
dashed histogram snows tne |stanway et al.| l |2007| ) distribution at 
z ~ 4 — 6. Lower: The resulting probability ot detecting TV lines 
with S/N> 10 in our sample, using the simulations described in 
the text. We observe no Lya with S/N> 10. 



2 out of 7 candidates. This estimate may be pessimistic, 
given the excellent photometric quality of the Hawk-I 
and WFC3 data, and the more careful cleaning of lower 
z interlopers compared to the V09 samples. However, 
the contaminant population may be changing at higher 
redshifts, and different and previously unstudied galaxy 
types may be entering the selection window. Ignoring 
these uncertainties, we repeated the Monte Carlo simu- 
lation for all possible choices of 5 candidates from our 7, 
finding that the probability of detecting no Lya line at 
S/N> 10 is still rather low, being typically 8%, and only 
in one case reaching 15% (this range depends on which 
candidates are excluded from the sample). 

Another explanation for the paucity of Lya detections 
would be physical evolution of either the galaxies or the 
surrounding IGM. The intrinsic strength of the Lya emis- 
sion is expected to increase at higher redshift as galax- 
ies become more metal- and dust-poor. The probability 
that these photons escape the galaxy and its surround- 
ings, however, depends on a series of complex (and not 
fully charact erized) phenomena in the IGM surrounding 
the galaxies (Zheng et al. 2010 and references therein), 
including the relative geometry and dynamics of gas and 
dust, e.g., backward scattering from wind-driven out- 
flowing shells (which can even boost the strength of the 
lines), or absorption by the damping wings of in- falling 
IGM along the line of sight. The presence of HI in prox- 
imity to the source can result in an absorption of the 
intrinsic Ly a by one order of magnitude or more ( Zheng 
et al. 2010), along with a broadening and redshifting ot 



the emerging line profile. Thus, one explanation for the 
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lack of Lya detections in our sample is a significant in- 
crease in the HI fraction of the IGM, xhi, at z ~ 7, lead- 
ing to a stronger absorption of the Lya flux. A similar 
effect could explain the obs erved decrease in the number 
of Lya emitters at z ~ 7 (Ota et al. 2010 Clement et 



al. , in prep. ), but results from these suryeys ar e still con 



tradictory flOuchi et al.||2010 
simulations ( |Dayal et a. "* 



1 



Hu et al ||2010| ). Detailed 
show that the IGM ab- 



sorption increases dramatically when the Universe is not 
fully ionized, leading to a significant absorption of the 
emerging Lya. The timescale of this effect around star 



formi ng galaxies is of the order of 100 Myr (Dayal et al. 
20101, shorter than the interval of cosmic time between 
z ~ 6. and z ~ 7. An additional prediction is that the 
asymmetry in the line profile is smoothed by the ve locity 
structure of the infalling IGM (|Dayal et aLl|2008[). Un- 



fortunately, the modest S/N in our only detection is too 
low to address this effect quantitatively. 

In conclusion, this work shows that the spectroscopic 
confirmation of z ~ 7 galaxy candidates is a challenging 
effort. However, these difficulties may not be only due to 
our current technological limits, but may also reflect the 
long-sought first evidence of the reionization process in 
the early Universe. Future surveys will definitely solve 
this fascinating puzzle. 



Observations were carried out using the Very Large 
Telescope at the ESO Paranal Observatory under Pro- 
gramme IDs 084. A-095, 181.A-0717. We thank an anony- 
mous referee for precious comments. We are grateful to 
P. M0ller and the whole ESO staff for their assistance 
during the execution of service observations. We ac- 
knowledge partial financial support by ASI. 



REFERENCES 



Bouwcns, R. J., Illingworth, G. D., Oesch, P. A., et al. 2010a, 
ApJ, 709, L133 

Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2010b, 
ApJ, 708, L69 

Bouwens, R. J., Thompson, R. I., Illingworth, G. D., et al. 2004, 
ApJ, 616, L79 

Bunker, A., Wilkins, S., Ellis, R., et al. 2009, ArXiv e-prints 
Castellano, M., Fontana, A., Boutsia, K., et al. 2010, A&A, 511, 
A20+ (C10) 

Dayal, P., Ferrara, A., & Gallerani, S. 2008, MNRAS, 389, 1683 
Dayal, P., Maselli, A., & Ferrara, A. 2010, ArXiv e-prints 
Douglas, L. S., Bremer, M. N., Lehnert, M. D., Stanway, E. R., & 

Milvang-Jensen, B. 2010, ArXiv e-prints (D10) 
Finkelstein, S. L., Papovich, C, Giavalisco, M., et al. 2009, ArXiv 

e-prints 

Giavalisco, M. 2002, ARA&A, 40, 579 

Giavalisco, M., Koratkar, A., & Calzetti, D. 1996, ApJ, 466, 831 
Greiner, J., Kriihler, T., Fynbo, J. P. U., et al. 2009, ApJ, 693, 
1610 

Hickey, S., Bunker, A., Jarvis, M. J., Chiu, K., & Bonfield, D. 

2010, MNRAS, 404, 212 
Hu, E. M., Cowie, L. L., Barger, A. J., et al. 2010, ArXiv e-prints 
lye, M., Ota, K., Kashikawa, N., et al. 2006, Nature, 443, 186 
Labbe, I., Gonzalez, V., Bouwens, R. J., et al. 2010, ApJ, 708, L26 
Mannucci, F., Buttery, H., Maiolino, R., Marconi, A., & Pozzetti, 

L. 2007, A&A, 461, 423 
McLure, R. J., Dunlop, J. S., Cirasuolo, M., et al. 2010, MNRAS, 

403, 960 

Oesch, P. A., Bouwens, R. J., Illingworth, G. D., et al. 2010, ApJ, 
709, L16 



Ota, K., Iye, M., Kashikawa, N., et al. 2010, ArXiv e-prints 
Ouchi, M., Mobasher, B., Shimasaku, K., et al. 2009, ApJ, 706, 
1136 

Ouchi, M., Shimasaku, K., Furusawa, H., et al. 2010, ArXiv 
e-prints 

Pentericci, L., Grazian, A., Scarlata, C, et al. 2010, A&A, 514, 
A64+ 

Reddy, N. A. & Steidel, C. C. 2009, ApJ, 692, 778 
Salvaterra, R., Delia Valle, M., Campana, S., et al. 2009, Nature, 
461, 1258 

Schaerer, D. & de Barros, S. 2010, ArXiv e-prints 

Shapley, A. E., Steidel, C. C, Pettini, M., & Adelberger, K. L. 

2003, ApJ, 588, 65 (S03) 
Stanway, E. R., Bunker, A. J., Glazebrook, K., et al. 2007, 

MNRAS, 376, 727 (S07) 
Stark, D. P., Ellis, R. S., Chiu, K., Ouchi, M., & Bunker, A. 2010, 

ArXiv e-prints 

Tanvir, N. R., Fox, D. B., Levan, A. J., et al. 2009, Nature, 461, 
1254 

Vanzella, E., Giavalisco, M., Dickinson, M., et al. 2009, ApJ, 695, 
1163 (V09) 

Wilkins, S. M., Bunker, A. J., Ellis, R. S., et al. 2010a, MNRAS, 
403, 938 

Wilkins, S. M., Bunker, A. J., Lorcnzoni, S., & Caruana, J. 

2010b, ArXiv e-prints 
Zheng, Z., Cen, R., Trac, H., & Miralda-Escude, J. 2010, ApJ, 

716, 574 



